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ABSTRACT 

We have performed 2D bulge/bar/disc decompositions using r and i-band images 
of a representative sample of nearly 1000 galaxies from the Sloan Digital Sky Survey. 
We show that the Petrosian concentration index is a better proxy for bulge-to-total 
ratio than the global Sersic index. We show that pseudo-bulges can be distinguished 
from classical bulges as outliers in the Kormendy relation. We provide the structural 
parameters and distributions of stellar masses of ellipticals, classical bulges, pseudo- 
bulges, discs and bars, and find that 32 per cent of the total stellar mass in massive 
galaxies in the local universe is contained in ellipticals, 36 per cent in discs, 25 per cent 
in classical bulges, 3 per cent in pseudo-bulges and 4 per cent in bars. Pseudo-bulges 
are currently undergoing intense star formation activity and populate the blue cloud 
of the colour-magnitude diagram. Most (though not all) classical bulges are quiescent 
and populate the red sequence of the diagram. Classical bulges follow a correlation 
between bulge Sersic index and bulge-to-total ratio, while pseudo-bulges do not. In 
addition, for a fixed bulge-to-total ratio, pseudo-bulges are less concentrated than 
classical bulges. Pseudo-bulges follow a mass-size relation similar to that followed by 
bars, and different from that followed by classical bulges. In the fundamental plane, 
pseudo-bulges occupy the same locus as discs. While these results point out different 
formation processes for classical and pseudo-bulges, we also find a significant overlap in 
their properties, indicating that the different processes might happen concomitantly. 
Finally, classical bulges and ellipticals follow offset mass-size relations, suggesting that 
high-mass bulges might not be simply high-mass ellipticals surrounded by discs. 

Key words: galaxies: bulges - galaxies: evolution ~ galaxies: formation ~ galaxies: 
fundamental parameters - galaxies: photometry - galaxies: structure 



1 INTRODUCTION 

The central components of disc galaxies have proven de- 
ceptively simple. The early view of galaxy bulges as scaled 
ellipticals is far from complete. Even the mere definition of a 
galaxy bulge is still prone to debate. Yet bulges hold crucial 
clues to galaxy formation and evolution. Current views dis- 
cern between bulges formed through violent processes, such 
as hierarchical clustering via minor mergers, named classical 
bulges, and those formed through longer time-scales, via disc 
instabilities and secular evolution processes, named pseudo- 
bulges (see IWvse et al.lll997l : iKormendv fc Kennicuttll2004 
for reviews). Given their dissimilar origins, it is naturally ex- 
pected that these different bulge categories should be struc- 
turally distinct. 

In fact, previous work have showed evidence that 
pseudo-bulges, when compared to classical bulges, tend 



E-mail: dimitri@mpa-garching.mpg.de 



to show younger stellar populations, kinematics supported 
by rotation, and less concentrated surface brightness pro- 
files, similar to those of discs, which can be parameterised 
by the corresponding Sersic indic e s (see e.g. CaroUo et al.l 
19971: iGadotti fc dos Anio3 l200ll: iKormendv et al.l l2006l : 



Drorv fc Fishej|2007l : iFisher fc Drorvl |2008| . and references 
therein). In relatively inclined galaxies, pseudo-bulges are 
seen to be more fiattened than classical bulges, and as flat- 
tened as discs. Some pseudo-bulges display a characteris- 
tic box/peanut shape if the galaxy is viewed edge-on. Such 
box/peanut bulges are believed to be the inner parts of bars 
tha t have evolved and b uckled off the plane of the disc (see 
e.g. lAthanassouial I2OO5I . and references therein). Such evi- 
dence is, however, limited. Since such studies brought along 
new perceptions, proofs of concept, and because pseudo- 
bulges are less conspicuous than classical bulges, they had to 
focus on a detailed analysis of small samples, often with high 
resolution data. For instance, pseudo-bulges often host nu- 
clear bars and/or nuclear spiral arms, which are difficult to 
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identify without images of liigli pliysical spatial resolution, 
or a careful analysis. 

Ellipticals would be formed through major mergers. The 
gas content in galaxies, as a dissipative component and the 
material to form new stars, is likely to play an important 
role in minor/major mergers and in secular evolution, and 
adds considerable complexity to these views. A more com- 
plete picture of galaxy formation and evolution demands 
a deep understanding of the consequences of such diverse 
processes, as well as the structural properties of the differ- 
ent galaxy components. It is still necessary to understand 
better how these different formation processes result in evo- 
lutionary and structural differences, how these differences 
translate into observables, and how this framework can be 
incorporated by theoretical developments. Questions such 
as how structurally different ellipticals, classical bulges and 
pseudo-bulges are, how different the host galaxies of classical 
and pseudo-bulges are, how often disc galaxies host pseudo- 
bulges as compared to classical bulges, how are their masses 
distributed, and in which conditions pseudo-bulges are im- 
portant, remain to be properly answered. Time is ripe for 
more extensive studies on the properties of bulges and ellip- 
ticals, taking advantage of large datasets, such as the Sloan 
Digital Sky Survey (SDSS), which are becoming more com- 
monly available. 

We have performed a detailed structural analysis of a 
sample of nearly 1000 galaxies from the SDSS. Through 
such a work, we were able to determine Sersic indices, scale- 
lengths and effective surface brightness for elliptical galax- 
ies, classical bulges and pseudo-bulges, and scale-lengths and 
central surface brightness for discs, taking into account the 
presence of a bar when necessary, as well as bulge-to-total 
and disc-to-total luminosity ratios, in three different band- 
passes {g, r and i). We have also determined colours and 
stellar masses for each galaxy component separately. Such 
structural parameters are explored throughout this paper, 
combined with several other physical properties obtained 
for the SDSS sample of galaxies. (A similar set of parame- 
ters were also obtained for bars, which are investigated in a 
sep arate paper.) 

lAlleii et al.l (|2006l ) and iBenson et all |20o3) have re- 
cently performed detailed structural analysis using large 
datasets. While their samples are substantially larger than 
ours, our methodology aims at a more careful analysis. For 
instance, we have include d a bar compon ent in the models 
fitted, since, as shown in iGadottil (|2008| ). if the bar is not 
taken into account in the modelling of barred galaxies, bulge 
properties can be significantly affected. We have produced 
fits in multiple bands for every galaxy, which have allowed 
us to estimate the colours of the different components. With 
such colours, we were able to directly estimate mass-to-light 
ratios for the different components, which results in more 
accurate stellar mass estimates. Finally, we have selected a 
sample of galaxies suitable for such image decompositions, 
and each fit was performed and checked individually. All 
this care assures us that the results from our analysis are 
trustworthy to carry the investigations we aimed for. 

A common practice in studies on bulges is to define 
pseudo-bulges as those with a Sersic index lower than a 
given threshold. However, the use of Sersic indices to dis- 
tinguish pseudo-bulges from classical bulges is likely to gen- 
erate considerable ambiguity, as the dynamical range of 



this parameter is small, compared to the uncertainty in 
individual measurem ents (see discussion in IGadottil l2008l : 
iDurbala et al."2008''). Inst ead, we have used the Kormendy 
relation "(jkormcridvl ll977l ) to identify pseudo-bulges. Some 
degree of structural similarity between classical bulges and 
ellipticals is borne out by such relation, which is followed 
by both kinds of systems. In contrast, pseudo-bulges, be- 
ing structurally different, do not follow a similar relation. 
Therefore, our identification of pseudo-bulges is not only 
more reliable but also better physically motivated. It is also 
important to point out that, due to the fact that our models 
include a bar component, the pseudo-bulges we identify in 
this study are not those pseudo-bulges associated with the 
box/peanut morphology. 

This paper is organised as follows. The next section de- 
scribes the sample selection and properties. Section [3] de- 
scribes the methodology implemented to produce the fits. 
Results are presented in Sect. |31 separated in several sub- 
sections. We discuss various aspects of our main findings in 
Sect.O including a comparison with previous work. Finally, 
we summarise our work and major conclusions in Sect. [51 
Throughout this work, we have used the following cosmo- 
logical parameters: Ho = 75 km s~^ Mpc~^, Qm = 0.3 and 
JIa = 0.7. Unless noted, photometric measurements refer to 
the i-band. All logarithms are base 10. 



2 SAMPLE SELECTION AND PROPERTIES 

We aimed for a sample which is concomitantly suitable for 
structural analysis based on image decomposition and a fair 
representation of the galaxy population in the local universe. 
We thus drew the sample from all objects spectroscopically 
classified as galaxies in the SDSS Data Release Two (DR2), 
applying the following criteria. We first selected all galaxies 
with redshift in the range 0.02 ^ z ^ 0.07. This provides a 
sample with a statistically meaningful number of galaxies, 
whose images have a relatively comparable physical spatial 
resolution. Images of closer galaxies have a level of detail 
much richer than those of galaxies at z ~ 0.05, as observed 
within SDSS. On the other hand, images of farther galaxies 
have a physical spatial resolution which is often not adequate 
to study their structural properties in some detail. 

Secondly, since dwarf galaxies are not object of our 
study, we excluded all galaxies with stellar masses be- 
low 10^° M(T). Galaxy s tellar masses were obtained from 
iKauffmann et al.l (l2003bh . At this stage, we have a volume- 
limited sample, i.e. a sample which includes all galaxies more 
massive than 10^° M© in the volume defined by our redshift 
cuts and the DR2 footprint. From now on, we will refer to 
this sample as our parent sample. 

Finally, we chose objects with an axial ratio b/a ^ 0.9, 
where a and b are, respectively, the semi-major and semi- 
minor axes of the galaxy, taken from the SDSS database 
(these are measured at the 25 mag arcsec"^ isophote in the 
p-band). This criterion assures us that our galaxies are very 
close to face-on, meaning that effects from dust attenuation, 
as well as projection effects, are minimised, and that bars 
can be distinctly seen, avoiding erroneous fits, and leading 
to a more reliable determination of the structural parame- 
ters and colours of the separate galaxy components. It also 
implies we do not need to worry about inclination correc- 
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Figure 1. Distributions of galaxy properties for the parent sample (top panels) and the final sample (bottom panels). From left to right: 
dust-corrected stellar mass, k-corrected 2-band absolute magnitude, D„4000 (corrected for [Nc in] contamination), ex tinction-corrected 

S O in] 5OO7A luminosity for AG N, and Petrosian concen tration index. The first three parameters were taken from iKauffmann et all 
2003d) . [O III] luminosities from IKauffmann et al.l ||20033), and concentration from th e SDSS database. The [O III] luminosity has been 
used in previous work as an indicator of the accretion rate onto the central black hole jHeckman et al.ll2004h . 



tions. These criteria resulted in a sample of 3375 objects. 
The images of each of these galaxies were individually in- 
spected to remove from the sample those galaxies which are 
either not truly face-on Q substantially disturbed by a com- 
panion or merger, overly faint or irregular, have images too 
close to the border of the CCD frame, as well as duplicate 
entries and those images where the presence of a bright star 
renders them useless. We also rejected galaxies with a < 4 
arcsec, as we deemed these galaxies too small for a detailed 
parametric image decomposition. This ensured that, for the 
vast majority of the sample, the effective radius of the bulge 
is larger than the PSF HWHM. Hence, reliable conclusions 
can be drawn even for small bulges (see discussion further 
below). Our final sample contains 963 galaxies, of which 407 
contain an active nucle us (AGN - mostly t ype 2) . AGN clas- 
sification is taken from IKauffmann et al.l 1 2003ah . 

Most of the galaxies in the final sample have 0.04 < 
z < 0.06. At this distance, 1 arcsec corresponds roughly to 
1 kpc. Since the SDSS spectra are taken through a circu- 
lar aperture of 3 arcsec in diameter, this means that, for 
the majority of our galaxies, the spectra obtained represent 
the central, bulge-dominated part of the galaxy, with little 
disc contamination. This will become clearer further below, 
where the effective radius of bulges are presented. Since the 
mean PSF FWHM in SDSS images is ^ 1.5 arcsec, the typ- 
ical physical spatial resolution in the images of our final 
sample is about 1.5 kpc. Most of these galaxies also have 
0.9 ^ b/a < 0.95. In many cases, this means that bulge and 
disc have slightly different geometrical properties, such as 
ellipticity e (i.e. 1 — b/a) and position angle. These differ- 
ences are helpful when fitting models to the images, since 
they provide further constraints to the models of the dif- 
ferent galaxy components. Finally, their typical stellar mass 
is similar to that of the Milky Way (i.e. « 5 x 10^° Mq - 
ISommer-Larsen fc Dolgovll2001^ . 

Let us now check whether our final sample is a fair rep- 

^ Axial ratios from SDSS measurements are not always correct. 
Galaxies with b/a clearly below 0.9 were removed from the sam- 
ple. 



resentation of the galaxy population in the parent, volume- 
limited sample, and whether our selection criteria to reach 
the final sample introduce any bias that have to be ac- 
counted for. With this aim, we plot in Fig. [1] the distri- 
butions of several relevant properties of the galaxies in both 
samples. One sees that these distributions are reasonably 
similar, in particular those of the AGN [O ill] luminosity. 
However, it is also apparent that the final sample includes 
a somewhat larger fraction of massive, brighter, quiescent 
and more concentrated galaxies. The more stringent crite- 
rion to reach the final sample from our parent sample is the 
axial ratio cut. Thus, we now investigate whether this crite- 
rion is the cause of such differences. The top panel of Fig. [2] 
shows the distributions of axial ratio for galaxies in the par- 
ent sample, separated according to their concentration C. 
From our structural analysis below, we find that most ellip- 
ticals, and few disc galaxies, have C > 3. Conversely, most 
disc-dominated galaxies, and few ellipticals, have C < 2.5. 
We thus use these thresholds to separate ellipticals and disc 
galaxies. It is evident that, due to their different intrinsic 
axial ratios, such galaxies have different distributions of the 
observed axial ratio, in such a way that the fraction of el- 
liptical galaxies increases at higher values of b/a. Therefore, 
it is indeed the axial ratio cut at b/a ^ 0.9 that produces 
the differences seen in Fig. [T] Some of the results presented 
below depend on the fraction of elliptical galaxies in the 
sample, and are thus affected by this selection effect. We 
can correct these results for such effect by comparing the 
fraction of galaxies with C > 3 in the parent sample, re- 
gardless of axial ratio, to that of galaxies with both C > 3 
and b/a J5 0.9. This provides us with a fairly good indication 
of how larger the fraction of elliptical galaxies in our final 
sample is compared to our parent, volume-limited sample. It 
turns out that this difference is a factor of 1.3. Such correc- 
tion factor will be applied when necessary with an explicit 
mention. 

Given our axial ratio cut, one might also question 
whether the results in this work concern galaxies which are 
intrinsically more eccentric than average. To check that, we 
have also selected 30 systems with C > 2.5 (i.e. mostly 
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Figure 2. Top: distribution of isophotal axial ratios for galaxies 
in the parent sample with concentration index below 2.5 (blue 
line) and above 3 (red line). Bottom: stellar mass plotted against 
redshift for all spectroscopically classified galaxies in SDSS DR2. 
Black lines indicate our cuts in mass and redshift. Blue and red 
dots correspond to galaxies with low and high dust-corrected z- 
band mass-to-light ratios, respectively. 



bulge-dominated and elliptical galaxies) and 0.5 < b/a < 
0.7, which otherwise comply with all other selection criteria 
of the final sample, and analyse them in a similar fashion. 
One could also ask if we are missing galaxies with very strong 
bars due to our cut in axial ratio. This would happen if the 
outer disc of the galaxy is not detected. We argue that this 
is not the case, based on the fact that the strongest bars 
we find, considering the bar-to-total luminosity ratio, are 
abo ut as strong as the s trongest ones found in other stud- 
ies (ISellwood fc Wilkinso n 1993: lReese et al.ll2007l : iGadottil 



l2008l : lDurbala et al.ll2008i ) 



Finally, we also need to check if there is a bias in our 
parent sample against galaxies close to both its low mass 
and high redshift cuts, which could occur due to magnitude 
limits in the survey. This is done by plotting galaxy mass 
against redshift for all spectroscopically classified galaxies 
in DR2 in the lower panel of Fig. [21 Since such a bias might 
depend on galaxy colour, we plot separately galaxies with 
low and high dust-corrected z-band mass-to- light ratio (from 
iKau ffmann et al. 2003b), using the peak of its distribution 
as boundary between the two groups. One sees that there is 
no clear need for such a volume correction, i.e. our parent 
sample includes effectively all galaxies more massive than 
10^° Mq in the defined volume. 

We have also compared the distributions of other phys- 
ical parameters, such as stellar surface density and metal- 
icity, r-band light-weighted age, number of neighbors, and 
the fluxes and equiv alent widths of important em ission lines 
in AGN, derived by IKauffmann et al.l l|2003al lbl |2004| ) and 



iGallazzi et al] l|2005h . in both our parent and final samples. 
These distributions are generally similar. We thus conclude 
that, after accounting for the effects introduced with our 
axial ratio cut, our final sample is a fair representation of 
the population of galaxies and AGN in the local universe 
with stellar masses larger than 10^^° Mq. This means that 
conclusions based on this sample are generally pertinent to 
galaxies that satisfy these conditions. 

Figure|3]shows examples of galaxies in our sample. Each 
row corresponds to a different galaxy category, concern- 
ing the central component. Elliptical galaxies are shown in 
the top row. The remaining rows show, from top to bot- 
tom, respectively, galaxies with classical bulges, galaxies 
with pseudo-bulges, and bulgeless galaxies. This grouping 
is based on the results from our parametric image decompo- 
sitions below. 



3 IMAGE DECOMPOSITION 

The fitting of each galaxy image in our sample was done 
using BUDDA v2.1, which is able to fit up to four different 
galactic components to a galaxy image, namely bulge, disc, 
bar, and a central source, usually adopted to account for the 
emission from bright, type 1 AGN. Details o n the code tech- 
nical a spec ts and capabilit ies can be found in lde Souza et all 
(|2004l ) and lGadottil (|2008l ). In particular, the latter discusses 
more specifically the treatment of the bar component. Here 
we briefiy recall the functions used to describe each compo- 
nent. 

The model compo nents are described by c oncentric, 
generalised ellipses (see I Athanassoula et al.lll990l ): 



(1) 



where x and y are the pixel coordinates of the ellipse points, 
a and b are the extent of its semi-major and semi-minor 
axes, respectively, and c is a shape parameter. The centre of 
the component, given by the coordinates xq and yo, is also 
fitted by the code, and is the same for all components fitted. 
Position angles and ellipticities were fitted by the code for 
every component, each one having a single position angle 
and ellipticity. When c — 2 one has a simple, pure ellipse. 
When c < 2 the ellipse is discy, and when c > 2 the ellipse 
is boxy. Bulges and discs were fitted using pure ellipses, but 
c was left free to fit bars, since bars are better described by 
boxy ellipses. 

Each component also has to follow a surface brightness 
radial profile. The disc profile is described by an exponential 
function (Freeman 1970), 



f^d{r) = MO + 1.086r-//i, 



(2) 



where r is the galactocentric radius, /lo is the disc central 
surface brightness, and h is the disc characteristic scale- 
length. 

The bu lge profile is describ ed by a Sersic function 
(|Sersidll968l . see lcaon et aLlll993l ), 



fibir) = fie + c„ 



(a 



l/n 



- 1 



(3) 



where is the effective radius of the bulge, i.e. the ra- 
dius that contains half of its light, fie is the bulge effec- 



Pseudo-bulges, classical bulges and ellipticals in SDSS 5 



• 


« 


• 


• 


• 












f 


# 






* 


■ m 


1 




■ 








# 

F 






• 

4 



Figure 3. Examples of galaxies in our sample. The four rows show from top to bottom, respectively: elliptical galaxies, galaxies with 
classical bulges, galaxies with pseudo-bulges, and bulgeless galaxies. For the three last rows, the first three columns from left show barred 
galaxies. This grouping is based on the results from our parametric image decompositions. 



tive surface brightness, i.e. the surface brightness at re, n 
is the Sersic index, defining the shape of the profile, and 
c„ = 2.5(0.868ri - 0.1 42). When n = 4 the profile is a de 
Vaucouleurs function (jde Vaucouleurd Il948l ). while n = 1 
corresponds to an exponential bulge. Elliptical galaxies were 
fitted with a single "bulge" component. 

The bar profile is also described by a Sersic function. 
For the bar, 



l/"bar 



(4) 



where Cn, bar = 2.5(0.868nbar — 0.142), and the other pa- 
rameters have similar definitions as for the bulge. Another 
bar parameter fitted by the code is the length of the bar 
semi-major axis, Lbar, after which the bar profile is simply 
truncated and drops to zero. 

Using the models obtained with budda, we have calcu- 
lated the absolute magnitude of each different galaxy com- 
ponent in all galaxies in the sample in the g, r and i bands. 
We have thus also calculated the integrated colours of each 
component separately. 

To avoid unreliable results that very often arise from 
automated image decomposition, we applied budda to each 
galaxy individually. Erroneous fits are particularly likely to 



arise, if not followed closely, when fitting barred galaxies. 
In such cases, the number of parameters fitted can be as 
high as 18, and the topology of the distribution of for 
the many possible models is very complex. We thus adopted 
a procedure to produce each fit, which is as follows. 

First, we needed to produce the image to be fitted. To 
this end, a cutout was created from the original SDSS galaxy 
image, with the galaxy in the centre, and containing enough 
sky background in order not to remove the galaxy outskirts. 
Usually, this means that the image size on a side is about 
twice the apparent diameter of the galaxy. Before fitting, 
each image was cleaned of spurious objects such as fore- 
ground stars. The sky background and PSF FWHM values 
given by SDSS were checked and corrected when necessary, 
which rarely occurred. When needed, the sky background 
was estimated by calculating the median pixel value of sev- 
eral areas free from objects in the galaxy image. When the 
PSF FWHM given by SDSS was deemed erron eous, we have 
made new estimates, by fitting a Moffat profile l|Moffatlll969l ) 
to a number of suitable stars near the galaxy, and choos- 
ing the median value found. The images were fitted without 
subtracting the mean background value. This is important, 
because a proper calculation of the x^y thus a reliable 
determination of the uncertainties, is based on the statistics 
concerning all photons that reached the detector. Evidently, 
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the value of the background is given to the code, and the 
model components are determined taking it into account. 

Second, before starting the fit, it is necessary to provide 
the code with initial guesses for all parameters to be fitted. 
The coordinates of the pixel with highest intensity were cho- 
sen as input to xo and yo. The geometrical parameters were 
roughly estimated by visual inspection of the image and the 
corresponding isophotal contours. When deemed necessary, 
these parameters were also calculated using the intensity 
moments in a few difi'erent circular apertures centred at the 
galaxy centre. Finally, the initial guesses for the parameters 
of the surface brightness profiles were estimated through the 
inspection of a pixel-by-pixel radial intensity profile of the 
galaxy. Tests with artificial images of unba rred disc galax- 
ies wi th high and low spatial resolution in Ide Souza et al.l 
showed that budda is generally able to reliably re- 
trieve the correct structural parameters. InXJadotti (2008), 
further tests with similar images showed that the uncertain- 
ties estimated by the code for each fitted parameter are real- 
istic. In Appendix [X] we perform new tests, and check how 
the fits in the present study depend on the input parame- 
ters, using SDSS r-band images of twenty barred galaxies in 
our sample. The choice for barred galaxies assures that the 
results from these tests concern typically difficult fits. The 
results are encouraging. They show that, in most cases, the 
code converges to similar results (generally consistent within 
the estimated uncertainties) for input parameters varying up 
to a factor of four. Furthermore, wrong fits can in general 
be recognised by their unusually high values. 

A critical point at this stage is to decide on which com- 
ponents the model of the galaxy should include. With this 
aim, the inspection of the galaxy image, isophotal contours 
and intensity profile was also very helpful. A bulge com- 
ponent was included in the model if the intensity profile is 
clearly not a single exponential function. In the absence of 
bars or spiral arms, it might not be clear only by inspecting 
the image if a galaxy has or not a disc component. The latter 
might become easier to identify when checking the isophotal 
contours for significant differences in position angle and/or 
ellipticity between the inner and outer parts of the galaxy. 
In addition, an outer exponential component in the inten- 
sity profile is naturally also an indication of the presence of a 
disc. A disc was included in the model whenever one of such 
signatures was found. To check for the presence of a bar, we 
first inspected the galaxy image. If in doubt, we then looked 
for the typical bar signatures in the isophotal contours and 
the intensity profile, i.e., respectively, an elongated structure 
with con stant position angle , and a fiat ledge in the profile 
(see e.g. iGadotti et"aLll2007^ . A bright AGN could be dis- 
cerned if the profile has a cusp in the centre. However, even 
bright AGN are smoothed out if the physi cal spatial reso - 
lution of the image is not high enough (see lGadottill200"^ . 
Accordingly, none of our images present such a feature, and 
thus no galaxy in our sample required an AGN component 
in the fit. 

Evidently, however, some cases are less clear, and one 
can not rule out the possibility of us fitting a truly ellipti- 
cal galaxy with a model containing bulge and disc, or con- 
versely, fitting an unbarred lenticular galaxy with only a 
bulge model. In some of the latter cases, a first fit resulted 
in a bad match to the galaxy image, usually a bulge with 
an obviously too low Sersic index. In most of these cases. 



another fit, including a disc component, produced a better 
match, revealing the true nature of the galaxy. To check that 
our decision on whether or not adding a disc component in 
a given fit is reliable, we have performed the following tests. 
We randomly selected 20 galaxies judged as ellipticals and 
produced test fits with a disc component. Conversely, we 
randomly selected 20 galaxies judged as unbarred lenticu- 
lars and produced test fits without a disc component. It 
turned out that, for the ellipticals, only one test fit resulted 
in a lower x^y and, for the lenticulars, two test fits resulted 
in similar values. Nevertheless, none of such cases are 
significantly improved fits, and most of these test fits re- 
sulted in unmistakably worse matches to the galaxy images. 
Furthermore, the disc component found in the tested el- 
lipticals never constributes with more than a few per cent 
of the total galaxy luminosity. It should also be mentioned 
that, due to the limited physical spatial resolution of our im- 
ages, we likely missed most bars shorter than Lbar « 2 — 3 
kpc, typically seen in very lat e-type spirals (later than Sc - 
lElmegreen fc ElmegreenI 19851 ). The whole of these faint bars 
is typically within 2 — 4 seeing elements, and thus they do 
not imprint clear signatures in either the isophotal contours 
or the intensity profile. Note that, in contrast, bulges with 
re of the order of one seeing HWHM can still be identified 
in the intensity profile, since they usually contain a much 
larger fraction of the galaxy light than these short and faint 
bars. 

This procedure was performed using several suitable 
tasks in IRAF0 The fits to each galaxy were first done in 
the r-band, until a satisfactory model is achieved. Then the 
same input parameters which resulted in the final r-band fit 
were used to fit the g and i-band images, with the luminos- 
ity parameters scaled accordingly. In addition, in some cases, 
geometrical parameters such as ellipticities and position an- 
gles were kept fixed at the values found with the r-band 
image. A model was deemed incorrect if any of the fitted 
parameters assumed a value which is clearly wrong, but this 
was uncommon. A typical case concerns the bulge Sersic in- 
dex. If n was found to be outside the range 0.8 ^ n ^ 6, then 
a new set of input parameters was tried in a new fit. Usually, 
the subsequent fits returned values for n within this range, 
as well as smaller uncertainties in the fitted parameters, and 
a lower x^- In these cases, we sticked to the new fit. In some 
other cases, however, the resulting model did not change af- 
ter changing the initial parameters. When this happened, 
we considered the model appropriate. Only occasionally a 
more detailed comparison between galaxy image and model 
was done, e.g. by comparing both intensity profiles. 

Figure |3] shows examples of the fits obtained for six 
galaxies in the sample. These include an elliptical galaxy, 
two galaxies with a classical bulge (one barred galaxy and 
one unbarred), two galaxies with a pseudo-bulge (again 
barred and unbarred), and a bulgeless galaxy. The distinc- 
tion between classical bulges and pseudo-bulges is explained 
in details in Sect. 14.21 One sees that the fits are generally 
quite good, even though the physical spatial resolution of 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 4. Examples of fits obtained. Eacli row refers to a single galaxy, showing, from left to right: the original i-band galaxy image, 
the model image obtained with BUDDA, the residual image, and surface brightness profiles. Horizontal lines in the galaxy images mark a 
length of 5 kpc at the galaxy distance. The first two images have intensities transformed to magnitudes per square arcsecond, and are 
displayed with the same levels in all rows. The residual images are obtained subtracting the model image from the original image, and 
are displayed with a much narrower range in intensity, which is also the same in all rows, in order to enhance residual sub-structures. In 
the residual images, darker pixels indicate where the galaxy is brighter than the model, whereas whiter pixels indicate where the model 
is brighter than the galaxy. The surface brightness profiles were obtained from cuts along the galaxy major axis, or along the bar major 
axis when there is a bar. The dashed line corresponds to the original galaxy image, whereas the black solid line corresponds to the model. 
When there is more than a single component, red, blue and green lines refer, respectively, to bulge, disc and bar. The dotted line in the 
bottom panel shows the residuals (galaxy — model). These results correspond to, from top to bottom, respectively: an elliptical galaxy, a 
galaxy with a classical bulge, a galaxy with a pseudo-bulge, a barred galaxy with a classical bulge, a barred galaxy with a pseudo-bulge, 
and a bulgeless galaxy. 



8 Dimitri A. Gadotti 



the images is restricted, and despite the fact that, given the 
size of our sample, we could not make a thourough com- 
parison between galaxy and model before accepting the fit 
to every single galaxy in the sample. This would require an 
unreasonable amount of time to complete the study for the 
whole sample. The discrepancies between galaxy and model 
are usually below 0.3 mag arcsec"^ and are mainly caused 
by sub-structures such as spiral arms, rings, lenses and bar 
ansae. It is clear that, if we can not ascertain that the fits are 
excellent for a single specific galaxy, the results from the fits 
are, on average, trustworthy, and thus suitable for statistical 
studies. 



4 RESULTS 

4.1 A comparison between different bulge-to-total 
indicators 

The morphological (Hubble) type of a galaxy is re- 
lated to several physical parameters that are an indica- 
tion of the galaxy form ation and evolution processes (e.g. 
[Roberts fc Haviieg|l994l ). One of the basic features on which 
the morphological classification of a galaxy hinges on is the 
relative importance of the bul ge component to the overall 
light distribution of the galaxy (|Hubblell 19261 ). Thus, it is of 
the utmost importance to be able to measure how relevant 
is the bulge component within a given galaxy in an objec- 
tive and quantitative fashion, even if there is some scatter in 
the relation between measured bul ge prominence and visu- 
ally determined Hubble type (e.g. iLaurikainen et al] l2007l : 
iGraham fc WorlevI l2008l ). Because it is difficult to objec- 
tively define the latter, the former is advantageous. 

In this work, since we have performed a detailed decom- 
position of each galaxy image into its main components, 
bulge prominence is readily accessible from the bulge-to- 
total luminosity ratio B/T, calculated from the models ob- 
tained with BUDDA. However, the employment of such tech- 
niques to very large samples of galaxies presents serious chal- 
lenges, as automated procedures are prone to large uncer- 
tainties. Thus, other parameters, which are not so difficult 
to determine, have been suggested as a measure of bulge 
prominence. One of such parameters is the Petrosian con- 
centration index, defined as the ratio i?90/-R50, where _R90 
and ii50 are, respectively, the radi i enclosing 90 and 50 per 
cent of the galaxy luminosity (see iKauffmann et al.ll2003bl . 
and references therein). Another possibility is to fit the en- 
tire galaxy as a single component, with a S ersic function, 
and th en use the Sersic index so obtained (e.g. lBlanton et al] 
l2005bl P: iHaussler et al] |2007| ) . 

Since we have measurements of B/T and bulge Sersic 
index for all galaxies in our sample, we can test how well 
concentration index and galaxy Sersic index relate to B/T, 
and how all these connected parameters relate to each 
other. This is done in Figs. [5] and |6l Concentration index 
is ta ken from the SPSS da tabase and galaxy Sersic index 
from iBlanton et al. I l|2005ah . The latter are updated val- 
ues after the findin g of a systematic bias in the values in 
iBlanton et alj (|2003l V Figure[5]shows that concentration cor- 
relates poorly with bulge Sersic index. Although one sees a 
tendency in the sense that bulges with higher values for 
n reside in galaxies with higher _R90/-R50, the spread is 
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Figure 5. Concentration index (R90/R50) plotted against bulge 
Sersic index n (top), bulge-to-total ratio B/T (middle) and the 
added contributions of bulge and bar to the total galaxy luminos- 
ity {B + Bar) /T (bottom). The correlation between concentration 
and n is remarkably poorer than that with B/T. Interestingly, the 
latter correlation is improved when the bar contribution is taken 
into account. Elliptical galaxies and bulgeless galaxies can be 
readily identified in the middle panel at B/T = 1 and B/T = 0, 
respectively. The red line in the bottom panel is a linear fit to 
the data points with < (B -I- Bar) /T ^ 0.6. For galaxies with 
B/T, or {B + Bar)/T, greater than fa 0.6, R90/R50 is roughly 
constant at Ri 3. 



very large. This might be partly due to the fact that n 
is related only to the bulge component, whereas i?90/-R50 
depends on the overall structure of the whole galaxy. An- 
other issue is that n can have high uncertainties compared 
to its dynamical range. Figure [5] also shows that _R90/-R50 
is well correlated with B/T. Interestingly, the correlation 
gets noticeably better when the bar contribution to the to- 
tal galaxy luminosity is taken into account. The correla- 
tion coefficient for _R90/i?50 x B/T is 0.82, while that for 
R90/R50x{B + Bar)/T is 0.86. Ahhough bars are more ex- 
tended structures than bulges, they are usually significantly 
less extended than discs, which explains the better correla- 
tion of R90/R50 with {B + Bar)/T. This means that there 
is an additional uncertainty when using i?90/-R50 as a proxy 
for B/T, as one galaxy with small bulge but massive bar can 
have a higher value of _R90/-R50 than an unbarred, but other- 
wise similar, galaxy. It is worth noticing that concentration 
rises linearly with {B + Bar)/T up to {B + Bar)/T « 0.6, 
after which it stays roughly constant at « 3, a value of con- 
centration similar to that of massive elliptical galaxies. A 
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Figure 6. Top left: galaxy Sersic index from iBlanton et al.l 
(|20053) plotted against our measurements of bulge Sersic index; 
top right: concentration index plotted against galaxy Sersic in- 
dex; bottom left: galaxy Sersic index plotted against B/T; bot- 
tom right: galaxy Sersic index plotted against {B + Bar)/T. The 
solid line in the top left panel depicts a perfect correspondence. 
The dashed line is a linear fit to all data points, fixing the slope 
of the fit to 1. Although there is a very large spread, it is clear 
that galaxy Sersic index is systematically higher than bulge Sersic 
index by, on average, ^ 0.8. 

linear fit to the data points with < (B -|- Bar)/T ^ 0.6 
gives 

i?90/i?50 = 1.93(±0.02) + 2.02(±0.05) x {B + Bar)/T. (5) 

Figure [6] shows how galaxy Sersic index relates to bulge 
Sersic index, concentration, B/T and {B + Bar)/T. All plots 
show a very large spread. In particular, galeixy Sersic index is 
systematically higher than bulge Sersic index by, on average, 
« 0.8. Moreover, galaxy Sersic index is a poor proxy for 
B/T. From the previous figure one sees that concentration 
index is much more reliable in this respect. The correlation 
coefficients for galaxy Sersic index as a function oi B/T and 
{B + Bar)/T are, respectively, 0.43 and 0.45. 

Nevertheless, it should be noted that both the galaxy 
and bulge Sersic indices in Figs.[S]and|S]are corrected for see- 
ing effects, whereas concentration index is not. It does seem, 
however, that see ing effects on concen tration index are small 
at low redshifts (|Blanton et al.l [20031 ). For a sample with a 
narrow range in redshift, as ours, such efi'ects seem to be 
practically negligible, since the systematic variation of the 
average physical size of the PSF is relatively small. We ver- 
ified this by plotting ii90/-R50 against redshift and finding 
no correlation whatsoever. Seeing effects would only intro- 
duce an uncertainty at each measurement, which is small 
compared to the spread in concentration seen at each B/T, 
or {B + Bar)/T, in Fig. [5] 

4.2 Identifying pseudo-bulges 

As mentioned in the Introduction, studies on the properties 
of pseudo-bulges are often characterised by having a physical 
spatial resolution higher than that of the SDSS images used 
here. It is thus legitimate to be concerned on how well can we 
in fact distinguish pseudo-bulges from classical bulges. One 
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Figure 7. Distribution of the ratio between bulge effective radius 
and PSF HWHM for all galaxies in the sample (excluding bulge- 
less galaxies). The two vertical dashed lines mark the positions 
where this ratio is 0.8 and 1. Only 3 per cent (10 per cent) of the 
bulges have effective radii below 0.8 (1) times the PSF HWHM. 

key aspect of this assessment is the comparison between the 
effective radius of the bulge rg and the PSF HWHM of the 
corresponding image. GadottJ l|2008h showed that the struc- 
tural properties of bulges can be reliably retrieved provided 
that re is larger than ~ 80 per cent of the PSF HWHM. This 
was gauged by fitting images of very nearby galaxies, at rel- 
atively high physical spatial resolution, and through the fit- 
ting of the same galaxies using artificially redshifted images. 
Figure [7] shows that only three per cent of our bulges do not 
satisfy this criterion. If one wants to be more stringent and 
extend this limit to 100 per cent of the seeing radius, one still 
sees that only 10 per cent of our bulges are below the thresh- 
old. We have not excluded these bulges from this study but 
have verified that their inclusion or removal does not result 
in any significant change in the results presented here. This 
gives us confidence that the structural properties we have 
obtained are generally reliable even for the less conspicuous 
bulges. Note that our decompositions are seeing-corrected. 

The top two panels in Fig. [8] show the mean effective 
surface brightness within the effective radius {/iejfl plotted 
against effective radius for the bulges and elliptical galaxies 
in our sample. Massive elliptical galaxies follow a well de- 
fined relation between (/ie) and r^, which is given by the solid 
line, constructed from a fit to the data. The two dashed lines 
are the fit's ±3a boundaries (over the zero-point and keep- 
ing the sl ope fixed). This r elation is known as the Kormendy 
relation (|Kormendvlll977h and is in fact one of the possible 
projections of the fundamental plane (we explore the funda- 
mental plane relations in more detail in Sect. 14.4] ). The slope 
of the Kormendy relation we find here is si milar (consis- 
t ent w ithin the uncertainties) to that found bv lPierini et al.l 
(|2003 ) using //-band images. Many bulges also follow this 
relation. If bulges are separated according to their Sersic in- 
dex n, one sees that most of the bulges with n J5 2 follow 
the ellipticals' relation, whereas most of those with n < 2 

^ Note that (/^e ) is not equal to /Xe , which is the effective surface 
brightness at the effective radius. In plots such as Fig. [8] and 
the fundamental plane (discussed further below), (/le) must be 
used, since it takes into account differences in the light (mass) 
distribution within different systems. 
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Figure 8. Graphical representation of our criterion to distinguish classical and pseudo-bulges. The two contiguous panels at the top 
show the mean effective surface brightness within the effective radius plotted against the logarithm of the effective radius. The two 
panels display elliptical galaxies and, separately, bulges with Sersic index above and below 2. Barred and unbarred galaxies are also 
indicated. The solid line is a fit to the elliptical galaxies and the two dashed lines point out its ±3(t boundaries. Bulges that lie below the 
lower dashed line are classified as pseudo-bulges, independently of their Sersic index. Bigger, white- filled circles represent systems with 
0.5 < b/a < 0.7, with similar colour coding. They do not show any clear distinct patterns. The separate panel at the top shows a density 
plot with all data points taken altogether. It shows that ellipticals, classical and pseudo-bulges indeed appear as independent groups in 
the Kormendy relation. The four histograms at the bottom show the distributions of Sersic index n and bulge-to-total luminosity ratio 
B/T for classical bulges and pseudo- bulges with n < 2 (left) and n ^ 2 (right). While the threshold in Sersic index can be considered as an 
approximation to identify pseudo-bulges, it is clear that it can generate many misclassifications. Our criterion can recognise more precisely 
pseudo-bulges, as those which are structurally different from classical bulges and ellipticals (also shown in the top right histogram). 



do not, lying below the relation. This is expected, since ex- 
ponential bulges arc often pseudo-bulges, which are believed 
to be structurally different from classical bulges and ellip- 
ticals (see e.g. ICarollclll999l : iKormendv fc Kennicuttll2004 
and references therein). However, many of the bulges with 
n < 2 do seem to follow the same relation of ellipticals. Con- 
versely, several bulges with n ^ 2 fall below this relation. It 
is thus clear that although the Sersic index threshold can be 



consider ed as an approximati on to identify pseudo-bulges 
(see e.g. lFisher fc Drorvlliooil ). it can introduce many mis- 
classifications. Since pseudo-bulges are expected to be struc- 
turally different from classical bulges, and to lie below the 
(/Je) —Te relation for ellipticals, we can use the lower dashed 
line in Fig.|S]to be the dividing line between classical bulges 
and pseudo-bulges. This means that pseudo-bulges satisfy 
the following inequality: 
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(Me) > 13.95 + 1.74 X log re 



(6) 



where (^e) and re are measured in the SDSS j-band, and 
is in parsecs. In fact, in this diagram, discs als o occupy the 
regio n below the relation defined by ellipticals (jPierini et al.l 
I2OO2I ). which is consistent with the expectation that pseudo- 
bulges are built from disc instabilities. From Fig. [8] one also 
sees that barred and unbarred galaxies behave in a similar 
fashion in what concerns their positions in the (fie) — re re- 
lation. The histograms in Fig. [8] show the distributions of n 
and B jT for classical bulges and pseudo- bulges with n < 2 
and n J5 2 in separate panels. They show that our criterion 
to separate these two bulge categories can identify bulges 
with n < 2, but with high values of B/T, as classical bulges. 
These bulges have often values of n close to 2. Conversely, 
we are also able to identify pseudo-bulges with n ^ 2. These 
bulges show often low values of B/T and also values for n 
close to 2. These results indicate that some bulges can be 
erroneously classified as classical or pseudo, if one uses the 
Sersic index for such classification, just for being by chance 
(due to measurement uncertainties) in the wrong side of the 
Sersic index threshold. One sees that, as expected, pseudo- 
bulges show a distribution of n unlike that of classical bulges 
and ellipticals. In addition, the distribution of n for classical 
bulges is wide, and ranges from values typical of ellipticals to 
values in the gap between pseudo-bulges and ellipticals. Fur- 
thermore, one also sees that classical bulges have a peak in 
B/T around 0.4, while the corresponding value for pseudo- 
bulges is about 0.1. 

Although it is clear that bulges with low and high Sersic 
index occupy diff'erent loci in the fie — re projection, Fig. [8] 
suggests that there is a more significant difference between 
all bulges and ellipticals than between the different bulges. 
To ch eck that, we have performed 2 D Kolm ogorov-Smirnov 
tests l|Fasano Sz Francc schini 198?! : IPress e t al. 1992l). The 
tests show that, in fact, in this projection, ellipticals dif- 
fer more from bulges (taken altogether) than bulges with 
high Sersic index differ from bulges with low Sersic index. 
In this context, we have checked, using the tests described 
in Sect.|3] that if a lenticular galaxy is mistaken by an ellip- 
tical, the shift in re could be at most only 0.04 dex. Thus, 
although such an effect, if present, would indeed bring el- 
lipticals closer to bulges in the Kormendy relation, its am- 
plitude is too small. In light of these results, one could ask 
whether a separation between classical and pseudo-bulges in 
the Kormendy relation is rather artificial. The density plot 
in Fig. [S] suggests that it is not: pseudo-bulges appear as 
an independent group of points, as do classical bulges and 
ellipticals. 

Previous studies also indicate that pseudo-bulges show 
star-forming activity, as opposed to classical bulges, which 
are constituted mainly by old stella r populations (e.g. 
iGadotti fc dos Aniod[200ll : iFishedl2006l ). We can check this 
difference withi n our sample using both the D„(4000) in- 
dex provided in lKauffmann et al.l (|2003bl ) and the SDSS fi- 
bre colour (i.e. the colour measured within the fibre aper- 
ture). The D„(4000) index is based on the 4000 A discon- 
tinuity seen in optical spectra of galaxies, and it is not 
sensitive to dust attenuation effects. Young stellar popula- 
tions have low values of D„(4000), as compared to old stel- 
lar populations. The fibre colours are integrated over the 
central region of the galaxy and can suffer from dust red- 
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Figure 9. Distributions of Dn(4000) (top) and fibre g — i colour 
(bottom) for pseudo-bulges and classical bulges. Tlie arrows in 
the bottom panel indicate an upper limit on how dust reddening 
could affect the colours of each bulge category, on average. 



dening. In Fig. (9] we show the distributions of D„(4000) 
and fibre g — i colour for pseudo-bulges and classical bulges. 
The result is striking. Our pseudo-bulges not only do in- 
deed show values of D„(4000) lower than those for classical 
bulges, indicating a higher level of star formation activity, 
but their distributions are also significantly apart (except 
for a tail of classical bulges showing star- forming activity). 
In fact, the peaks of the D„(4000) distributions for pseudo- 
bulges and classical bulges coincide with the peaks in the 
bimodal distribution of D„(4000) seen in both our parent 
and final samples (see F ig. [1}. Such bimodality was also 
found bv lAUen et al] (|2006l . who concluded that it is pro- 
duced by bulge-dominated galaxies on the one hand, and 
disc-domin ated galaxies on the other hand. More recently, 
iDrorv fc^ ishcr (2007) suggested that such bimodality is a 
result of the two different bulge families: pseudo-bulges are 
hosted by galaxies with high overall star formation activity, 
whereas classical bulges reside in galaxies that underwent vi- 
olent relaxation processes in the past, such as mergers, and 
thus have formed most of their stars in these past events. 
Both suggestions are not mutually exclusive, since galaxies 
with pseudo-bulges are often disc-dominated, as opposed to 
bulge-dominated galaxies, which often host classical bulges. 
In fact, our results indicate that the most disc-dominated 
galaxies almost always host pseudo-bulges, which are almost 
always under current intense star-forming activity, whereas 
classical bulges are mostly quiescent and reside mainly in 
more bulge-dominated galaxies. Interestingly enough, star- 
forming classical bulges form a group of points of their own 



12 Dimitri A. Gadotti 



in the density plot of Fig. [S] between the lower dashed line 
and the solid line. We will come back to these points further 
below. We have verified that if one separates bulges using 
their Sersic index with a threshold at n = 2 the resulting 
D„(4000) distributions are similar, and one also sees a tail 
of star-forming bulges with n J5 2. In this case, however, the 
peaks of such distributions are slightly closer, and the dis- 
tribution for bulges with n < 2 is broader than that of our 
pseudo-bulges. This gives support to our criterion to identify 
pseudo-bulges. 

The colour distributions are also separated, though in 
a less dramatic way. Pseudo-bulges show bluer colours than 
classical bulges, on average, but the difference in the mean 
colour values amounts only to 0.2 mag. Evidently, dust red- 
dening does play a role in making these distributions more 
similar, as compared to those of D,i(4000). Some dust atten- 
uation is expected even in face-on galaxies (e.g. Driver ct al. 
I2OO8I ). The fact that the colour distribution of pseudo- bulges 
is wider than that of classical bulges is also a result of dust 
attenuation. On the one hand, there are those pseudo-bulges 
with not much dust attenuation, whose intense star forma- 
tion shows up as the resulting blue colours (and the blue 
tail in the distribution). And, on the other hand, there are 
those pseudo-bulges with dust-obscured star formation that 
produce the red colours seen in the red tail of the distribu- 
tion. One can have an estimate of the dust reddening in our 
b ulges, using estimates of t he total galaxy dust attenuation 
in lKauffmann et al.l (j2003bl ) as an upper limit (i.e. consider- 
ing that all dust attenuation in the galaxy strikes the bulge 
light). We calculated the average dust att enuation in the i- 
band u sing the results for the 2-band in iKauffmann et aP 
(|2003bl ). and assuming that dust attenuation goes as A~°'^ 
(where A is the observed wavelength) . The average i-band at- 
tenuation we thus find for our elliptical galaxies and galaxies 
with classical and pseudo-bulges is, respectively, 0.04, 0.16 
and 0.48 mag. This corresponds to g — i colour excesses of, 
respectively, 0.02, 0.08 and 0.15 mag. The colour excess for 
classical and pseudo-bulges are indicated by the arrows in 
the lower panel of Fig. [9] Interestingly, the difference in the 
position of the peaks of the D„(4000) distributions amounts 
to 0.60, while the corresponding difference for fibre g i 
colour is, as mentioned above, 0.20 mag (without any dust 
correction) and 0.27 (with the upper limit dust correction). 

Since the D„(4000) values and fibre colours are taken 
from SDSS spectra, which are themselves taken through a 
fixed&ove aperture of 3 arcsec (thus 1.5 arcsec in radius), and 
since pseudo-bulges are often smaller than classical bulges, 
one might ask whether the striking difference in D,i(4000) 
between pseudo-bulges and classical bulges in our sample is 
only due to aperture effects. One could argue that, due to 
the smaller extent of pseudo-bulges within the fibre aper- 
ture, the D„(4000) values (and similarly the fibre colours) 
for pseudo-bulges are contaminated with star formation in 
the disc. In fact, the average effective radius of pseudo-bulges 
and classical bulges is, respectively, 0.82 arcsec and 1.14 arc- 
sec. However, one must bear in mind that bulge light dom- 
inates over disc light through ~ 2 times the bulge effective 
radius from the galaxy centre. This is borne out in our sur- 
face brig htness profiles, and h a s also been found b y others 
(see e.g. iMoriondo et al.lll998l : iMorelh et alll2008l 'l. Never- 
theless, to be reassured that the difference seen in Dn(4000) 
is not just a spurious effect, we have verified that the differ- 
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Figure 10. Comparison between structural parameters of ellipti- 
cal galaxies and bulges obtained in the g and i-bands. Top: Sersic 
index; bottom: effective radius. Solid lines depict a one-to-one cor- 
respondence. Bulges are separated according to our classification 
between classical and pseudo-bulges. The absence of systematic 
effects indicates that dust does not play a significant role in the 
overall results of the decompositions. 



ence in the mean values of D„(4000) for pseudo-bulges and 
classical bulges is similar when all bulges are separated in 
bins of similar effective radius. In particular, we have verified 
that all results from Fig. hold if one considers only bulges 
with 7-e > 1 arcsec, or Ve > 1.5 arcsec, where little disc con- 
tamination is expected. Furthermore, we have also verified 
that the SDSS fibre colours are in fact much more similar to 
the integrated bulge colours, as calculated with the deter- 
mined BUDDA models for the g and i-band images, than the 
corresponding disc colours. We thus conclude that one can 
not blame disc contamination within the SDSS spectra to 
explain the higher star formation activity and bluer colours 
seen in our pseudo-bulges. 

Our sample was designed in a way to minimise dust 
effects. In fact, most of our galaxies have dust attenuation 
measurements in the low end of the distribution of dust at- 
te nuation for the w h ole sam ple of more than 10^ galaxies 
in IKauffmann et aL I (|2003bl ). Nevertheless, it is clear that 
the high star formation activity seen in pseudo-bulges is ac- 
companied to some extent with dust obscuration. Could this 
somehow affect the structural parameters obtained? To an- 
swer this question we plot in Fig. [10] bulge structural pa- 
rameters obtained with the g-band images against the same 
parameters obtained with the i-band images. Dust has more 
pronounced effects in bluer bands and thus, if there is some 
systematic effects caused by dust in our decompositions, 
these should show up in these plots, in particular for pseudo- 
bulges. For instance, dust attenuation could in principle re- 
sult in systematically lower values for n, as well as smaller 
values for r^, in the g-hand, as compared to those in the i- 
band. Figure [10] shows no substantial systematic effects. We 
thus conclude that our decompositions are not significantly 
disturbed by dust. T his is consistent with the findings in 
iMoriondo et al.l (|l998l ) , who found that bulge structural pa- 
rameters do not change significantly with galaxy inclination. 
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Figure 11. Bulge Sersic index plotted against bulge- to- total lu- 
minosity ratio (top left), bulge effective surface brightness (top 
centre) and bulge effective radius (top right). And bulge-to-total 
luminosity ratio plotted against bulge effective surface brightness 
(bottom left) and effective radius (bottom right). Blue, green and 
red circles represent, respectively, pseudo-bulges, classical bulges 
and ellipticals, with b/a ^ 0.9. Bigger, white-filled circles with 
similar colour coding represent systems with 0.5 < b/a < 0.7. 
They follow similar patterns. 

4.3 Scaling relations 

Figure [TTI shows the relations between bulge Sersic index and 
bulge-to-total luminosity ratio, effective surface briglitness 
and effective radius for pseudo-bulges, classical bulges and 
ellipticals. The relations between bulge-to-total luminosity 
ratio and effective surface briglitness and effective radius 
are also shown. It is striking how these three stellar systems 
populate well-defined and separate regions in some of these 
diagrams. Systems with larger n tend to be more conspic- 
uous and extended. One sees a correlation between n and 
B /T for classical bulges but n ot as much for ps e udo-b ulges, 
corroborating the results from lFisher fc Drorvl (|2008l ) . who 
arrived at a similar conclusion using a smaller sample. It is 
worth noticing that n rises with re for bulges but it is rather 
constant for ellipticals. Furthermore, there is a correlation 
between B/T and Ve for both classical and pseudo-bulges, al- 
though they seem to follow offset relations. Ellipticals tend 
to be more extended than bulges and have lower surface 
brightness than classical bulges. Note that no difference is 
seen in Fig. [TT] when comparing galaxies with different ax- 
ial ratios. We have also found that systems with larger n 
and B/T tend to be in more massive galaxies, although the 
scatter is large. 

In Fig.[l2]we show that, as expected, galaxies with more 
extended bulges or discs tend to be more luminous and mas- 
sive. This is also true for systems with lower axial ratios. 
One might also ask if a similar relation is found if one con- 
siders the component mass, rather than the total mass in 
the galaxy. Although several stellar systems seem to follow 
such a relati on, there might be exceptions, such as globular 
clusters (see iBurstein et al.l[l997l : iKissler-Patig et aLllioO^ : 
iBarmbv et al.ll2007l ). This is answered in Fig. 1131 One sees 
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Figure 12. Disc scale-length and bulge effective radius plot- 
ted against galaxy 2-band absolute magnitud e (k-corrected to 
z = 0.1 ) and dust-corrected stellar mass, from iKauffmann et al.l 
ll2003d) . Blue, green and red circles represent, respectively, galax- 
ies with pseudo-bulges, galaxies with classical bulges and ellipti- 
cals, with b/a ^ 0.9. Bigger, white-filled circles with similar colour 
coding represent systems with 0.5 < b/a < 0.7. They follow sim- 
ilar patterns. 

that such a correlation indeed exits for pseudo-bulges, clas- 
sical bulges, bars and discs, separately, regardless of galaxy 
axial ratio. This figure also shows fits to the data, where 
scale- lengths go oc A/", where a is a constant that defines 
the slope of the relation and M is the component stellar 
mass. The values we have found for a for bars, discs, pseudo- 
bulges, classical bulges and ellipticals are, respectively, 0.21, 
0.33, 0.20, 0.30 and 0.38, with a mean la uncertainty of 0.02; 
thus a seems to increase for more massive components. It is 
interesting to note that bars and pseudo-bulges share a very 
similar value of a (and these are the only components that 
do so), which is different from that of classical bulges at a 
5(T level. This not only suggests a different formation mech- 
anism for classical and pseudo-bulges, but also reveal a close 
connection between pseudo-bulges and bars. Such a connec- 
tion is expected if pseudo-bulges indeed form through disc 
instabilities. Nonetheless, despite such difference in slope, 
the locus occupied by classical bulges in such diagram con- 
tinues rather smoothly from the locus of pseudo-bulges. This 
suggests that although there might be two different dom- 
inant formation processes for classical and pseudo-bulges 
there can be a significant number of cases where the differ- 
ent processes are at play concomitantly. Thus, the sequence 
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Figure 13. Scale-lengths of elliptical galaxies, classical bulges, 
pseudo-bulges, discs and bars, plotted against the correspond- 
ing component stellar mass, as indicated. Solid lines are fits 
to the data. Bigger, white-filled circles represent systems with 
0.5 < b/a < 0.7, with similar colour coding. They show a similar 
behaviour. 

drawn by pseudo-bulges and classical bulges could repre- 
sent a sequence where the dominant process changes rather 
smoothly from one end to the other. 

Another striking feature in Fig. 1131 concerns ellipticals. 
They are more offset from classical bulges than classical 
bulges are offset from pseudo-bulges, as 2D Kolmogorov- 
Smirnov tests indicate. In fact, one sees that most classical 
bulges can not be considered as ellipticals that happen to be 
surrounded by a disc, at least not ellipticals more massive 
than 10^*^ M0. How robust this result is against systematic 
effects from the decompositions? We explore now two possi- 
bilities. The first is that some of the ellipticals could actually 
be lenticulars whose discs were not identified in the prepa- 
ration of the fits. As mentioned in Sect.[3l we produced test 
fits in which a disc component was forced in 20 ellipticals, 
and found that in only one case the inclusion of a disc could 
be justified. However, such a disc contributes with only a 
few per cent of the total galaxy luminosity, and the result- 
ing shift in re, towards lower values, is only of the order of 
a few hundred parsecs. The second possibility is that the 
input values given to the code to fit are systematically 
larger for ellipticals. We have thus performed yet another 
test, again producing test fits for 20 ellipticals but this time 
with the input values for re cut by half. In most cases, the 
result was similar to the original fit (which incidently at- 
tests the robustness of the code). In a few cases, the shift 
in the final value for is again only a few hundred parsecs 
towards lower values, and in all such cases the resulting 
is substantially worse than that of the original fit. We thus 
conclude that such systematic effects could in principle bring 
some of the low-mass ellipticals closer to classical bulges in 
Fig. 1131 but the result that at least a significant fraction 
of classical bulges can not be considered as (massive) ellip- 
ticals surrounded by discs seems to be robust. Note that 
iFisher fc Drorvl (|2008l ) argue that classical bulges resemble 
low-luminosity ellipticals. 

With these data, we can test previous claims that 
the Hubble sequence is scale-free, i.e. that the relation be- 
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Figure 14. Ratio between the scale-lengths of bulges and discs 
Te/h plotted against B/T for pseudo-bulges and classical bulges. 
Although there is some scatter, it is evident that both quantities 
arc correlated, which rules out the concept of a scale-free Hubble 
sequence. Interestingly, pseudo-bulges and classical bulges seem 
to follow offset relations. For the same B/T, pseudo-bulges have 
on average higher values of Ve/h. Bigger, white-filled circles with 
similar colour coding are pseudo-bulges and classical bulges in 
galaxies with 0.5 < b/a < 0.7, and portray similar patterns. 

tween the scale-len gths of bulges and di scs does not vary 
with Hubble type dCourteau et al.lll996l ). In this context, 
iGraham fc Prietol (|l999[ l showed that fixing the bulge Sersic 
index, when fitting galaxy surface brightness profiles, can 
lead to an overestimate or an underestimate of re, depend- 
ing on the true value of n. Since n is a free parameter in our 
decompositions, our results are not subject to such system- 
atic biases. In Fie. 1141 we plot Ve/h against B/T for pseudo- 
bulges and classical bulges. This figure shows clearly that 
Ve/h and B/T are correlated, and thus the Hubble sequence 
is not scale-free, if it is a B/T sequence. iBalcells et al] (|2007D 
and lLaurikainen et al.l l|2007l ) reached similar results, though 
with much smaller samples. Figure [14] also indicates that the 
relations for pseudo-bulges and classical bulges are offset: for 
a constant B/T, pseudo-bulges have on average higher val- 
ues of re//i. 

It is well known that galaxies can be separated depend- 
ing on the locus they occupy in a colour-magnitude diagram. 
Some galaxies follow a tight colour-magnitude relation, in 
the sense that more luminous galaxies are redder, known 
as the red sequence, whereas others spread beside the red 
sequence, orig inating a s warm of points known as the blue 
cloud (see e.g. lGladderTe t al. 1998; Bell ct al. 2007;, and ref- 
erences therein). Such a diagram is a powerful tool, able to 
provide clues on the formation and evolu ti onary processes 
of gal axies (see e.g. iDe Lucia et al.ll2007^ ■ iDrorv fc Fisher! 
(2007'), based on a careful analysis of a sample of 39 galax- 
ies, suggested that galaxies with pseudo-bulges populate the 
blue cloud, while galaxies with classical bulges populate the 
red sequence. But colour can be subject to dust attenua- 
tion effects, and galaxy magnitude is actually used as a 
proxy to galaxy stellar mass. Since we have measurements of 
D„(4000), an efficient, dust-insensitive index to identify stel- 
lar populations with different mean ages, as well as carefully 
modelled stellar masses, we can go a step further and plot 
D„(4000) against galaxy stellar mass, as is done in Fig. [151 It 
should be noted that often one uses the total galaxy colour in 
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Figure 15. The 4000A break index plotted against galaxy dust- 
corrected stellar mass. Virtually all elliptical galaxies are found 
in the red sequence, while practically all galaxies with pseudo- 
bulges, and bulgeless galaxies, are found in the blue cloud. Galax- 
ies with classical bulges can be found in both the red sequence 
and the blue cloud, although the majority of them occupy the red 
sequence. 



the colour-magnitude diagram, whereas our D„ (4000) values 
correspond to the central galaxy colour. That is, however, 
an asset, since, as argued in Sect. 14.21 the light that goes 
through the SDSS fibre, from which Dn(4000) is measured, is 
in our sample mostly emitted by stars in the bulge, meaning 
that our Fig. [15] is mainly a comparison between ellipticals, 
classical bulges and pseudo-bulges, avoiding in most cases ef- 
fects resulting from disc light. Nonetheless, we have verified 
that similar results emerge if we plot the colour-magnitude 
diagram in its usual form, using dereddened g — i Petrosian 
colou rs from the SDSS databa se and absolute magnitudes 
from iKauflmann et al.1 (|2003bl ). In such case, however, the 
blue cloud is not as detached from the red sequence as in 
Fig. 1151 most likely due to dust attenuation. Our results 
show clearly that virtually all elliptical galaxies are found in 
the red sequence, and practically all galaxies with pseudo- 
bulges, and bulgeless galaxies, are found in the blue cloud. 
However, although most of the galaxies with classical bulges 
do indeed populate the red sequence, a significant fraction 
of them are found in the blue cloud (see also Fig. [9} . 



4.4 The fundamental plane 

Elliptical galaxies follow well-defined relations between their 
effective radius re, mean effective surface brightness (/ie) and 
central velocity dispersion ao, known as the fundamental 
plane (FP). Furthermore, it is possible to show that such re- 
lations arise in virialisod systems (Djorgovski & Davis 1987; 
lDressler et~al.ll 19871 : iBender et al.llT992l '). Thus, the FP pro- 
vides im portant clues on the formation processes of stellar 
svstems. iBender et al.l (|l992l ) proposed a formulation of the 
FP, defining a space, whose three axes, ki, k,2 and ac3, are 
defined in terms of Ve, (fie) and ao, and are proportional 
to important physical parameters: ki is proportional to the 
logarithm of the mass, K2 is proportional mainly to the log- 
arithm of the surface brightness, and ks is proportional to 
the logarithm of the mass-to-light ratio. Here, mass means 



dynamical mass, i.e. stellar mass plus the dark matter mass 
content. 

What is the locus occupied by pseudo-bulges and clas- 
sical bulges in the FP? We can answer this question using 
Te and (/ie) from our decompositions, and velocity disper- 
sion measurements from the SDSS database. We have cho- 
sen to use the estimates provided in the Data Release 6 
(DR6), which are based on results from a direct fit to the 
spectra. In earlier releases, the estimates are an average of 
results from direct-fitting and Fourier-fitting, but the lat- 
ter systematically overestimates the velo city dispersion for 
low-mass galaxies (|Bernardi et al. I l2003al . their Fig. 17). In 
fact, we have verified that the velocity dispersion estimates 
from DR2 are systematically larger (by ~ 10 — 20 per cent) 
than those from DR6, for the galaxies in our sample, when 
the velocity dispersion is below ~ 150 km/s. Due to the 
limited spectral resolution of the SDSS spectra, it is recom- 
mended to use only spectra with signal-to-noise ratio above 
10, and without warning flags. We have verified that only 
one galaxy in our sample does not comply with the former 
criterion, and only 8 do not comply with the latter. Before 
we can use these measurements we have to apply aperture 
corrections due to the fact that these velocity dispersions are 
obtained through a fixed fibre aperture of 3 arcsec, which 
is physically different from galaxy to galaxy, depending on 
their dist ances. We th us used the prescription provided by 
[jorgensen et al.l l|l995l ) and converted the SDSS velocity dis- 
persion measurements to the velocity dispersion at 1/8 of r^, 
(jg. This empirical prescription is based on measurements for 
elliptical and lenticular galaxies, and it is unclear if it is also 
valid for disc-dominated galaxies. However, the corrections 
applied are small, and do not affect our results significantly. 

We thus calculated ki, K2 and K3 for elliptical galax- 
ies, cla ssical bulges a nd p seudo-bulges, using the defini- 
tions in lBender et al.l (|l992l ). and applying an offset in sur- 
face b right ness to facilitat e com parison with iBender et al.l 
(| 19921 ) and iBurstein et al.1 (|l997l i. who used measurements 
in the -B-band (the same procedure was done also by 
iBernardi et al.ll2003bl) . The top panels of Fig. [TS] show that 
the relation we have found for against ki , c onsidering only 
elliptica l galaxies, is similar to that found bv IBernardi et akl 
(|2003bl ). which used SDS S i-band data for a samp le of nearly 
9000 early-type galaxies. IBernardi etU] (|2003bh used esti- 
mates of the velocity dispersion that include Fourier-fitting 
results, which, as mentioned above, overestimate the velocity 
dispersion for low-mass galaxies. Accordingly, their relation 
is somewhat less steep than ours, even though they have 
excluded estimates below 100 km/s. We find that 



= 0.27fi:i - 0.11. 



(7) 



In addition, classical bulges also seem to follow the same 
relation, although somewhat offset to higher K3 values, re- 
gardless of the galaxy bein g barred or unbarred (see also 
iFalcon-Barroso et al.1 I2OO2I '). C lassical bulges f ollow m ore 
closely the relation found by IBernardi et al.l (|2003bl ). It 
should be noted that their sample also includes early-type 
disc galaxies. To select their sample of early-type galaxies 
they have included only galaxies with concentration index 
-R90/-R50 larger than 2.5. This does not mean, however, 
that their sample contains only elliptical galaxies. In fact, 
we have checked that 99 per cent of our ellipticals have 
-R90/-R50 > 2.5, but this is also the case for 76 per cent of our 
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Figure 16. Elliptical galaxies, classical bulges and pseudo-bulges in the K-space formulation of the fundamental plane (FP). The top 
three panels show the edge-on view of the FP, K3 plotted against fti , while the bottom panel shows its nearly face-on view, K2 plotted 
against k\ . Barred and unbarred galaxies are also indicated. In the top panels, the solid line is a fit to our ellipticals, while the dashed 
line is the fit obtained by Bernardi et al ] l|2003bl ') for nearly 9000 early-type galaxies. The dashed line in the bottom panel shows the 
limit of the zone of avoidance. Arrows indicate how some important physical parameters vary accross the FP. Bigger, white-filled circles 
represent systems with 0.5 < b/a < 0.7, with similar colour coding. 



galaxies with classic al bulges, and 8 per cen t of our galaxies 
with pseudo-bulges. iBernardi et al.l (l2003al ) also used other 
criteria to avoid including late-type galaxies in their sam- 
ple, but it is unlikely that these criteria excluded most of 
the galaxies with classical bulges. This can also explain why 
our relation for ellipticals is slightly different from theirs. 

Pseudo-bulges in both barred and unbarred galaxies lie 
significantly and systematically above these relations. This 
offset can not be attributed to effects due to dust attenua- 
tion alone. Even considering the upper limit for dust atten- 



uation discussed in Sect. 14.21 such a correction will bring 
the points down vertically by 0.19 in K3, while one sees 
that the offset is 0.3, on average. In this projection, 2D 
Kolmogorov-Smirnov tests indicate that pseudo-bulges dif- 
fer more from classical bulges than ellipticals differ from 
classical bulges or bulges altogether. Pseudo-bulges also he 
in a locus in the x ki projection of the FP which is differ- 
ent from those of both ellipticals and classical bulges (bot- 
tom panel of Fig. I16|l . Again, correcting for upper limit dust 
effects would put pseudo-bulges up along the k,2 axis only 
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Figure 17. Distributions of fiQ, h and D/T for all disc galaxies in the sample. The first two panels also show the values of the median 
and standard deviation of each distribution, as well as the mean Itr uncertainty in individual measurements. Bin sizes are fa 2a. 



by 0.27. In this projection, the statistical tests indicate that 
pseudo-bulges differ from classical bulges as much as ellipti- 
cals differ from classical bulges or bulges altogether. In both 
projections, pseudo-bulges lie in region s close to those oc - 
cupied by disc-dominated galaxies (see lPierini et al.ll2002l ). 
We note that barred galaxies tend to be more separated 
from their unbarred counterparts in the case of galaxies with 
pseudo-bulges, as compared to galaxies with classical bulges. 
Barred galeixies seem to have slightly higher values of both 
K2 and K3, on average. This is also confirmed with the statis- 
tical tests. Finally, it is also worth noting that, as expected 
from the results in Sect. 14.31 systems with lower axial ratios 
seem to generally share the same properties as those with 
b/a ^ 0.9. 

4.5 Discs 

Figure fTTl shows the distributions of the disc parameters (/iq, 
h and D/T). It is worth noticing that the distribution of 
disc-to-total ratio shows two peaks which match the peaks 
in the bulge-to-total ratio distributions for classical bulges 
and pseudo-bulges (shown in Fig. |8} . The two peaks in the 
distribution of D/T are at about 0.5 and 0.9. Bars account 
typically for 10 per cent of the galaxy luminosity. In Fig. [TS] 
we show the scaling relation between disc scale-length and 
central surface brightness. When one separates galaxies by 
their bulge type (classical vs. pseudo), one sees that, on aver- 
age, galaxies with pseudo-bulges have more extended discs, 
with fainter central surface brightness, as compared with 
galaxies hosting classical bulges, albeit with significant over- 
lap. A Kolmogorov-Smirnov test rejects the null hypothesis 
of no difference between such distributions at a 92 per cent 
confidence level. It is also interesting to ask how the stellar 
mass in discs varies between these two galaxy classes. The 
stellar mass in discs of galaxies with pseudo-bulges is, on av- 
erage, 1.15 X lO^'' M0, whereas that in discs of galaxies with 
classical bulges is 1.93 x lO^'' M©, thus a factor of 1.7 larger. 
A similar factor (1.8) corresponds to bars, when present. We 
can also ask how do their bulge masses compare. The stel- 
lar mass within pseudo-bulges is typically 2.2 x lO'' M©, 
whereas in classical bulges this is 1.41 X 10^" Mq (see Ta- 
ble 14)1, thus a factor of 6.4 larger. Hence, in terms of stel- 
lar mass, there is more similarity between discs (and bars) 
in galaxies with pseudo-bulges and those in galaxies with 
classical bulges, than between classical and pseudo-bulges 
themselves. Furthermore, the mean difference between the 
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Figure 18. Scaling relation between disc scale-length and cen- 
tral surface brightness. The histograms at the bottom show 
that galaxies with pseudo-bulges have more extended discs, with 
fainter central surface brightness, on average, as compared with 
galaxies hosting classical bulges, albeit with significant overlap. 
A Kolmogorov-Smirnov test rejects the null hypothesis of no dif- 
ference between such distributions at a 92 per cent confidence 
level. 

total galaxy mass of galaxies hosting pseudo-bulges and of 
those with classical bulges is a factor of 2. Therefore, if one 
scales down a typical galaxy with a classical bulge, dividing 
the mass of each of its components by 2, one ends up with a 
galaxy that has a disc (and a bar) with mass similar to that 
of the corresponding component in galaxies with pseudo- 
bulges. However, the bulge of such scaled-down galaxy is 
around 3 times more massive than the typical pseudo-bulge. 
If we assume that the physical processes that lead to the 
formation of discs and bars are the same regardless of the 
resulting bulge, then this difference suggests that the for- 
mation processes of classical bulges and pseudo-bulges are 
typically distinct. 

4.6 A compendium of galaxy structural 
parameters 

The full results from our decompositions are given in Ta- 
bles [1] and (2] Table [3] shows the parameters describing the 
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Table 1. Structural parameters of bulges, discs and bars in our sample from j-band image fits. Galaxy identifications, uncertainties, 
and similar tables in the g and r bands are also available. The printed version of the paper contains only a sample. The full table is 
available in the electronic issue of the journal. 
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0.046 


20.17 


1.84 


20, 


,94 


2 


,02 


4, 


,18 


20, 


,56 


1.73 


0.63 


0.74 


5.08 


2, 


,4 


0.58 


0.29 


0.12 


1 


.54 


0, 


,61 


0.042 


20.78 


3.25 


21, 


,36 


0, 


,58 


1 


,17 





,0 


0.0 


0.0 


0.0 


0.0 


0, 


,0 


0.03 


0.96 


0.0 





.91 


0, 


,71 


0.049 


0.0 


0.0 


20, 


16 


1 


,48 


3, 


,17 





,0 


0.0 


0.0 


0.0 


0.0 


0, 


,0 


1.0 


0.0 


0.0 


1 


.53 


0, 


,55 


0.047 


20.80 


3.91 


21, 


,07 


2 


11 


3, 


,01 


20, 


,01 


1.52 


0.62 


0.86 


4.15 


2.36 


0.0 


0.0 


0.0 


1 


.76 


0, 


,95 


0.046 


20.33 


2.6 


21. 


,66 


1 


,77 


3, 


,88 





,0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.34 


0.65 


0.0 


2 


.02 


0, 


,94 


0.045 



Columns (1) and (2) show the disc central surface brightness and scale-length, respectively. Columns (3), (4) and (5) give, respectively, bulge effective 
surface brightness, effective radius and Sersic index. Bar parameters are given in columns (6) to (11), respectively: effective surface brightness, effective 
radius, ellipticity, Sersic index, semi-major axis and boxyness. Bulge-to-total, disc-to- total and bar-to-total luminosity ratios are given in columns (12), 
(13) and (14), respectively. The value of the fit is shown in column (15). Finally, columns (16) and (17), show, respectively, the seeing HWHM of 
the observation and the galaxy redshift. Intensities arc in units of mag arcsec""^ and lengths in kpc. Seeing is in arcsec. 



Table 2. Absolute magnitudes and integrated colours of bulges, discs and bars in our sample from the budda models. The printed 
version of the paper contains only a sample. The full table is available in the electronic issue of the journal. 
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0.0 


0.0 


1.22 


0.38 


1 


.00 
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0.0 


0.0 
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-21.34 


-19, 


,40 


-20, 


,68 
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.26 


0.0 


0.0 


0.0 
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0.46 


1 


.86 


0.58 


0.0 


0.0 


-20.17 


-20.64 


-20.93 


-15 


,81 


-16, 


,86 


-17, 


.23 


0.0 


0.0 


0.0 


0.76 


0.29 


1 


.42 


0.37 


0.0 


0.0 


0.0 


0.0 


0.0 


-20, 


,67 


-21 


,39 


-21 


.73 


0.0 
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0.0 


0.0 


0.0 


1 


.06 
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0.0 


0.0 


0.0 


0.0 


0.0 


-20, 
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0.0 


0.0 
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,17 
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,13 


-20, 
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0.0 


0.0 


0.0 
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1 
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0.0 
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0.0 


-19, 


,86 


-20, 


,79 


-21 


.33 


0.0 


0.0 


0.0 


0.0 


0.0 


1 


.46 
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0.0 


0.0 


0.0 


0.0 
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,15 
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0.0 
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.44 
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-16, 


,75 
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,25 
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0.0 
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-19, 


,26 
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.12 
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,51 
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0.0 
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1 


.24 
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0.0 


0.0 


0.0 


-18, 


,99 


-19, 


.80 


-20, 
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0.0 


0.0 


0.0 


0.0 


0.0 


1 


.15 


0.35 


0.0 


0.0 


0.0 


0.0 


0.0 


-19, 


,66 


-20, 


,31 


-20, 


,70 


0.0 


0.0 


0.0 


0.0 


0.0 


1 


.03 
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0.0 


0.0 
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-19.80 


-19, 


,46 


-20, 


.22 


-20, 


,53 


-17.72 


-18.74 


-18.82 


0.96 


0.32 


1 


.07 


0.30 


1.09 


0.07 


-19.67 


-20.16 


-20.39 


-15, 


,23 


-16, 


,32 


-16, 


,74 


0.0 


0.0 


0.0 


0.72 


0.23 


1 


.51 


0.42 


0.0 


0.0 


0.0 


0.0 


0.0 


-19, 


,17 


-20, 
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,40 


0.0 
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0.0 


0.0 


1 


.23 
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0.0 


0.0 
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-19, 
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-20, 


,35 


-20, 


,62 


-17.87 
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-19.14 


1.27 


0.49 


1 


.28 


0.26 


1.27 


0.32 


-19.39 


-20.19 


-20.34 


-19, 


,01 


-17, 


,96 


-19, 


,66 


0.0 


0.0 


0.0 


0.94 
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.65 


1.70 


0.0 


0.0 



Columns (1) to (3) show the disc absolute magnitude in the 9, r and i bands, respectively. Columns (4) to (9) show the same parameters for bulge 
and bar, respectively, as indicated. Columns (10) and (11), (12) and (13), and (14) and (15) give the integrated (9 — i) and (r — i) colours of disc, 
bulge and bar, respectively, as indicated. 



distributions (median and standard deviation), as well as 
the mean la uncertainty in individual measurements, of /ie, 
Te and n for elliptical galaxies, classical bulges and pseudo- 
bulges. 

Evidently, theoretical work on galaxy formation and 
evolution should comply with these results, either by pre- 
dicting compatible structural properties, or by using these 
results as a starting point to build realistic models. The 
same thing can be said about the masses of the different 
components, and their mass distributions, which are pre- 
sented in the next section. These results are also useful to 



test automated procedures of bulge/disc/bar 2D decompo- 
sition, with no checking or direct human intervention at the 
decomposition of every single object. An automated code 
able to reproduce these results using the same input data 
(i.e. the same images) should be deemed successful and reli- 
able to be used with much larger samples, which are lately 
becoming commonly available. The results from our image 
decompositions are available for download|f| These include 

■* Seelhttp://www. mpa-garching.mpg.de/~dimitri/buddaonsdss/| 
buddaonsdss.html . 
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Table 3. Parameters describing the distributions (median and standard deviation), as well as the mean 
1(T uncertainty in individual measurements, of /Je , re and n for elliptical galaxies, classical bulges and 
pseudo-bulges from the i-band decompositions. 





median 


/Xe (mag arcsec ^) 
sd 


Icr 


median 


re (kpc) 
sd 


Icr 


median 


n 
sd 


1(7 


Elliptical 


20.8 


0.4 


0.10 


3.0 


1.0 


0.15 


3.8 


0.9 


0.5 


Classical 


19.6 


0.4 


0.05 


1.0 


0.4 


0.03 


3.4 


1.3 


0.4 


Pseudo 


20.5 


0.5 


0.07 


0.7 


0.3 


0.03 


1.5 


0.9 


0.2 



Table 4. Parameters describing the distribu- 
tions (median and standard deviation) of the stel- 
lar mass in elliptical galaxies, classical bulges, 
pseudo-bulges, discs and bars. 





median (10^° Mq) 


sd (10"' Mo) 


Ellipticals 


5.15 


4.37 


Classical 


1.41 


1.98 


Pseudo 


0.22 


0.99 


Discs 


1.66 


1.62 


Bars 


0.40 


0.61 



all structural parameters of bulges, discs and bars in the 
g, r and i bands, absolute magnitudes and integrated g — i 
and r — i colours of each component separately, as well as 
model images for the whole galaxy and for each component 
separately. 

4.7 The stellar mass and luminosity budgets at 

2 ~ 

With the g — i integrated colour of each galaxy compo- 
nent, we have obtained the corresponding mass-to-light ratio 
in the i-band, using the results found by iKauffmann et al.l 
(|2007l . see their Fig. 15). They showed that the mass-to- 
light ratio in the i-band is well correlated with g — i colour 
over a range of more than a factor 10 in mass-to-light ratio 
and with a scatter of less than 0.1 dex. They found that 
the g ~ i colour is the colour that shows t he tightest rela- 
tion w ith mass-to-light ratio. Furthermore, iBell fc de Jond 
l|200ll ) showed that the combined effects of stellar age and 
dust attenuation are such that the stellar masses obtained 
here with these mass-to-light ratios are not significantly af- 
fected by dust. 

With the mass-to-light ratio and the total luminosity of 
each component we were thus able to calculate the stellar 
mass content in each component separately. Interestingly, we 
have verified that the total mass in each galaxy, calculated 
by adding the masses of the individual galaxy components, 
agrees v ery well with the galaxy stellar mass estimated in 
[Kaufl^ mann et al.l (|2003bl '). with a scatter of less than 0.02 
dex. The parameters that characterise the stellar mass dis- 
tributions of each galaxy component (median and standard 
deviation) are shown in Table |4] All such distributions are 
skewed towards lower mass values. 

We were also able to compute the fraction of the total 
stellar mass content in our galaxies contained within each 
component separately, thus obtaining the stellar mass bud- 
get in the local universe for massive galaxies (more massive 
than 10^" Mq). The following fi gures take into account the 



selection efi'ect due to our axial ratio cut discussed in Sect. 
[5] Elliptical galaxies contain 32 per cent of the total stel- 
lar mass, and the corresponding values for discs, bulges and 
bars are, respectively, 36, 28 and 4 per cent. Classical bulges 
contain 25 per cent of the total stellar mass, while pseudo- 
bulges contain 3 per cent. Hence, 11 per cent of the stellar 
mass in bulges are in pseudo-bulges, the remaining being 
in classical bulges. The uncertainty in these fractions from 
Poisson statistics only is between 1 and 2 percentage points. 

The corresponding luminosity fractions are as follows: 
27 per cent of the total luminosity from the galaxies in our 
sample (in the i-band) come from elliptical galaxies, 46 per 
cent from discs, 24 per cent from bulges, and 3 per cent from 
bars. Classical bulges emit 20 per cent of the total luminos- 
ity, while pseudo-bulges account for 4 per cent. Hence, 16 
per cent of the luminosity from bulges come from pseudo- 
bulges, the remaining coming from classical bulges. 

As expected, luminosity ratios, such as the bulge-to- 
total ratio B/T, change in different bands. For pseudo- 
bulges, the median value of B/T in the (;-band is 0.053, 
whereas in the i-band it is 0.075. For classical bulges, the 
median value of B/T in the p-band is 0.349, and in the i- 
band 0.377. The median value of Bar/T also increases from 
0.087 in the p-band to 0.095 in the i-band. The median value 
of D/T decreases from 0.702 in the g-band to 0.684 in the 
i-band. 

It is also worth noting the fraction in number corre- 
sponding to the different galaxy categories. From the total 
number of galaxies in our sample, 22 per cent are ellipticals, 
76 per cent are disc galaxies with bulges, and 2 per cent are 
bulgeless galaxies (i.e. pure discs). Within disc galaxies with 
bulges, 68 per cent host classical bulges, while 32 per cent 
have pseudo-bulges. Bars are found in 42 per cent of all disc 
galaxies. 



5 DISCUSSION 

5.1 Comparison with previous work 

IXii en et al.l l|2006t ) also performed parametric decomposition 
of a large sample of galaxy images. As mentioned in the In- 
troduction, their sample is significantly larger than ours, 
but our methodology aims at a more careful analysis. They 
found that the stellar mass content in the local universe is 
distributed as 13 per cent in elliptical galaxies, 58 per cent in 
discs, 26 per cent in classical bulges, 1.5 per cent in pseudo- 
bulges and the rema inder in low-luminosity blu e spheroids 
(|Driver et al.|[2007al . see also lBenson et al]|2007l who found 
that 35 — 51 per cent of the stellar mass is in discs). The 
figures do not change substantially when dust is taken into 
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account jPriver et al.ll2007bl ). A direct comparison with our 
results is difficult, since their sample was not selected with 
a cut in stellar mass. As we avoided low- mass galaxies, it is 
natural to expect that we find a larger fraction of the local 
stellar mass in ellipticals, and a corresponding lower frac- 
tion in discs. Interestingly, however, their fraction of stellar 
mass in classical bulges is similar to what we find. This is 
somewhat unexpected, considering the different sample se- 
lection criteria, since their sample should contain a larger 
fraction of disc-dominated galaxies than our sample. This 
might be explained by the a bsence of a bar component in 
their models. As discussed in iGadottil (j2008l l. it is possible 
to obtain good (but wrong) fits for barred galaxies, with 
no bar component in the model, because the bulge model 
absorbs the bar. Thus, their bulge-to-total luminosity ra- 
tios could be overestimated. They have considered pseudo- 
bulges as those with B/T < 0.5 and Ve/h < 0.5, and found 
them to be generally bluer than the remaining bulges. From 
our Fig. 1141 one sees that such criteria still include many 
classical bulges. However, they found a stellar content in 
pseudo-bulges smaller than ours, albeit within the uncer- 
tainties. A possible explanation for that could be that a 
number of their pseudo-bulges were smoothed out by the 
PSF. Although they have secured that the galaxies in their 
sample have effective radii l arger than the PSF HWHM, it 
was shown in IGadottil |2003) that it is the effective radius of 
the bulge that has to be considered , and not that of the whole 
galaxy. This issue is aggravated in iBenson et all (|2007l ). who 
used SDSS images up to a redshift of 0.3. 

In or der to try a more m eaningful comparison with the 
results in lDriver iiaDiioozi) we can assume that all galax- 
ies with stellar masses below 10^^° Mq are pure discs. This 
certainly underestimates the stellar mass content in bars 
and pseudo-bulges, but can give us a reasonable assessment 
of the distribution of stellar mass in discs, ellipticals and 
classical bulges, as if we have m ade no mass cut. Using the 
stellar mass function derived bv lBell et al.l (|2003h . using g- 
band data for a local sample of SDSS galaxies, we find that 
22 per cent of the stellar mass in the local universe is in 
galaxies with stellar masses below 10^° Mq. If aU this mass 
is in discs, then our local stellar mass budget becomes as 
follows. Discs contain 50 per cent of the stellar mass, ellipti- 
cals contain 25 per cent of stellar mass, classical bulges and 
pseudo-bulges contribute with, respectively, 20 per cent and 
2 per cent, and finally bars contain 3 per cent of the stellar 
mass. One thus sees a more reasonable agreement, although 
we find a somewhat larger mass fraction in ellipticals. 

In IGadottil { 20081 ) it was argued, with indirect means, 
that the stellar mass content in bars at z ^ is » 12 pe r 
cent of the total stellar mass (see also IWeinzirl et al]|2008h . 
which is a factor of three larger than what we find here. This 
discrepancy can be explained, at least partially, by two fac- 
tors. Firstly, the bar fraction assumed in the previous work 
is 70 per cent, which is typic ally found in studies dedicated 
to estimate this value (see e.g. lEskridge et ai]|2000l ). whereas 
the fraction of disc galaxies with bars in this work is 42 per 
cent. Secondly, we have selected our galaxy sample avoiding 
galaxies with low stellar mass. Thus, the stellar mass budget 
we find concerns massive galaxies, and it is thus biased in 
favour of ellipticals, which can not contain the bars we are 
addressing here. The fact that our bar fraction is low can 
also be explained by two factors. First, as discussed in Sect. 



(3] we miss most of the short bars, due to the limited spatial 
resolution in SDSS images. Second, bar fraction seems to be 
larger in low-m ass galaxies than in more massive galaxies 
l|Barazza et al.li200a V 

5.2 Classical bulges with star formation 

One interesting finding from Figs. |5] and [TJ] is that while 
pseudo-bulges show almost always intense star formation, 
classical bulges can either be quiescent or forming young 
stars. This indicates that galaxies with classical bulges can 
be further divided according to their values of D„(4000). 
We have thus separated classical bulges in two different cat- 
egories, namely those with D,i(4000) < 1.7, and those with 
D„(4000) ^ 1.7, and find that these categories have typically 
different bulge-to-total ratios and stellar masses. Galaxies 
hosting classical bulges with D„(4000) < 1.7 have median 
values oi B/T = 0.26 and mass equal to 0.8 x 10^° Mq, 
whereas those with D„(4000) ^ 1.7 have median values of 
B/T = 0.42 and mass equal to 1.8 x 10^" Mq. This sug- 
gests that the former are dominated by early-type spirals 
whereas the latter are mostly lenticular galaxies. To verify 
that, we have inspected randomly selected samples of these 
two galaxy categories and concluded that, in fact, galax- 
ies hosting classical bulges with D„(4000) < 1.7 show more 
often spiral structure and blue discs than galaxies hosting 
classical bulges with D„(4000) ^ 1.7 (see Fig. [T9j. In addi- 
tion, we show in Fig. [20] that these galaxies have markedly 
different distributions of the concentration index _R90/i?50. 
Furthermore, galaxies with pseudo-bulges also show a dis- 
tribution of _R90/-R50 distinct from that of both categories 
of galaxies with classical bulges. Elliptical galaxies, however, 
have a i?90/i?50 distribution more similar to that of galax- 
ies with classical bulges and D„(4000) ^ 1.7. We thus see 
three peaks in the di stribution of con c entra tion. A related 
finding is reported in iBailin fc Harri^ (|2008l ). As expected, 
the median Sersic index of star-forming classical bulges is 
2.8, thus lower than that of quiescent classical bulges, for 
which it is 3.6. 

A plausible explanation for the existence of galaxies 
hosting classical bulges with conspicuous star formation ac- 
tivity is the coexistence of a classical bulge and a pseudo- 
bulge in the same galaxy. This would naturally occur if, 
after the formation of the classical bulge (e.g. via a mi- 
nor merger), the galaxy disc develops an instability (such 
as a bar) able to induce the observed central star formation 
through processes of secular evolutionlfl Since pseudo-bulges 
are about an order of magnitude less massive than classical 
bulges (see Table 14)1, such composite bulges would have the 
structural signature of a classical bulge, concurrent with the 
star formation activity of a pseudo-bulge, as observed. The 
fact that the D„(4000) values of star- forming classical bulges 
and pseudo-bulges are similar (see Fig. [5} is consistent with 
this picture. Interestingly, we find that 34 per cent of galax- 
ies with classical bulges have D„(4000) < 1.7 (50 per cent of 

^ Another possibility is that star-forming classical bulges are 
those caught during a minor merger that enhances central star 
formation and will eventually help building the bulge. However, 
note that we have rejected galaxies with disturbed morphologies 
(Sect. [I. 
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Figure 19. Examples of galaxies hosting classical bulges with D„(4000) ^ 1.7 (top) and Dn(4000) < 1.7 (bottom). It is evident that 
the latter show more often spiral structure and blue discs. 
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Figure 20. Distribution of concentration (R90 / R50) for pseudo- 
bulges, classical bulges with D,i(4000) < 1.7, classical bulges with 
Dn(4000) ^ 1.7, and elliptical galaxies. 

which hosting bars), the remainder having D„(4000) ^1.7 
(with a fraction of barred galaxies of 40 per cent). Such 
excess of bars among star-forming classical bulges is also 
generally consistent with the picture of a composite bulge. 
In this context, it is worth noting that classical bulges with 
star formation lie between pseudo-bulges and classical bulges 
without star formation in the density plot of Fig. [8] (i.e. be- 
tween the lower dashed line and the solid line) , and also have 
intermediate masses, as mentioned above (see also Table |4]). 
They thus bridge pseudo-bulges and classical bulges without 
star formation in the mass-size relation in Fig. 1131 If such 
a picture is correct, then the Hubble classification could be 
understood as mainly a sequence that goes from elliptical 
galaxies successively to disc galaxies with classical bulges, 
composite bulges and pseudo-bulges. 

5.3 Formation of bulges and ellipticals 

The results in Sect. |3] all corroborate previous ideas that 
bulges, as identified through their photometrical properties, 
are present in two main classes, with distinct structural 



properties, indicating different formation processes. Those 
more similar to elliptical galaxies, namely classical bulges, 
seem to be formed through violent processes, such as hier- 
archical clustering, whereas those dissimilar bulges, namely 
pseudo-bulges, seem to be still under formation processes, 
more likely as a result of their host galaxies' slow internal 
evolution, which involves disc instabilities and angular mo- 
mentum transfer. 

These inferences are borne out by the fact that pseudo- 
bulges, when compared to their classical counterparts, have 
currently intense star formation activity (Fie:s.l^and ll5|l . are 
less concentrated at fixed bulge-to-total ratio (Fig. I14p . and 
occupy the same locus in the FP as discs (Fig. I16|) . Bulge 
Sersic index correlates with bulge-to-total ratio for classi- 
cal bulges but not for pseudo- bulges (Fig. Illf) . In addition, 
pseudo-bulges follow a mass-size relation similar to that fol- 
lowed by bars, and different from that followed by classical 
bulges (Fig.llSp. Furthermore, the discs (and bars) in galax- 
ies that host pseudo-bulges are not as different from those in 
galaxies that host classical bulges as pseudo-bulges are dif- 
ferent from classical bulges, in what concerns their relative 
stellar masses (Sect. 14. 5p . This suggests that disc and bar 
formation are similar in both galaxy classes, whereas bulge 
formation is distinct. 

However, we have also found indications that the pro- 
cesses that lead to the formation of classical and pseudo- 
bulges might happen concomitantly. In this context, the 
processes that lead to the formation of a pseudo-bulge are 
dominant at one end, while those that lead to the formation 
of a classical bulge are dominant at the other end, across 
the bulge population. This is suggested by the significant 
overlap in the distributions of bulges with n > 2 and n < 2 
in the Kormendy relation (Fig. |8]), which is in fact statisti- 
cally more significant than the overlap in the distributions 
of ellipticals and bulges altogether, as discussed in Sect. 14.21 
We have also found classical bulges with young stellar pop- 
ulations, typical of pseudo-bulges (Fig. [9]). Furthermore, al- 
though classical and pseudo-bulges have significantly differ- 
ent slopes for their mass-size relations, such relations join 
smoothly, with pseudo-bulges at the low-mass end, and clas- 
sical bulges at the high-mass end (Fig. I13|l . 

Classical bulges seem to bear some structural similarity 
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to ellipticals, since they follow the same Kormendy relation 
(Fig. |8} , and occupy a similar locus in the edge-on view of 
the FP (Fig. I16|l . However, they also follow offset mass-size 
relations (Fig. ll3p . This could in principle challenge the idea 
that classical bulges are simply ellipticals that happen to be 
surrounded by a disc. It should be noted, though, that this 
result concerns only massive ellipticals. Low-luminosity el- 
lipticals could, if surrounded by a disc, resemble classical 
bulges, but this can not be evaluated with the present data. 
We have verified that accounting for systematic effects could 
bring some of the low-mass ellipticals closer to the relation 
followed by classical bulges. Our tests indicate that such 
systematic effects should be restricted to about 10 per cent 
of the ellipticals. In the worst case scenario that all ellipti- 
cals are affected by systematic effects, the offset between the 
mass-size relations of classical bulges and ellipticals would 
still be present, but restricted to high- mass systems. Thus, 
the result that high-mass bulges can not be considered as 
high-mass ellipticals surrounded by discs seems to be ro- 
bust. We note that in the face-on view of the FP fFig. I16|) . 
pseudo-bulges, classical bulges and ellipticals occupy three 
significantly distinct loci, which should provide clues to the 
different formation and evolution processes. 



gests that the different formation processes might happen 
concomitantly, with different processes being dominant in 
different cases. 

We found that classical bulges and ellipticals follow off- 
set mass-size relations, suggesting that high-mass bulges can 
not be considered as high-mass ellipticals that happen to be 
surrounded by a disc. 

We provided distributions and typical values (includ- 
ing uncertainties in individual measurements) of all struc- 
tural parameters obtained (with the exception of bar pa- 
rameters, which will be addressed in a separate paper). We 
calculated the stellar mass content and distribution for each 
galaxy component, and showed that, considering galaxies 
more massive than 10^'' M©, in the local universe, 32 per 
cent of the total stellar mass is contained in ellipticals, and 
the corresponding values for discs, bulges and bars are, re- 
spectively, 36, 28 and 4 per cent. Classical bulges contain 
25 per cent of the total stellar mass, while pseudo-bulges 
contain 3 per cent. We also provided the corresponding lu- 
minosity and number fractions. In particular, we find that 
approximately a third of disc galaxies hosts pseudo-bulges. 
Such figures should be reproduced by successful models of 
galaxy formation and evolution. 



6 SUMMARY AND CONCLUSIONS 

We have determined several structural parameters for ellip- 
tical galaxies, bulges, discs and bars through reliable multi- 
band image fitting of a representative sample of nearly 1000 
local, massive galaxies in the SDSS. We showed that the 
Petrosian concentration index is a better proxy for bulge- 
to-total ratio than the global Sersic index. 

We showed that, while bulge Sersic index can be consid- 
ered as a criterion to distinguish pseudo-bulges from classical 
bulges, a more reliable, and physically motivated, separation 
can be made using the Kormendy relation. While classical 
bulges follow the relation set by ellipticals, pseudo-bulges 
do not, independently of their Sersic index. This shows that 
classical bulges bear some structural similarity to elliptical 
galaxies, while pseudo-bulges are structurally different. 

Using D„(4000), which is insensitive to dust attenua- 
tion, we demonstrated that pseudo-bulges are currently un- 
dergoing intense star formation activity, and that virtually 
all pseudo-bulges are located in the blue cloud of the colour- 
magnitude diagram. In contrast, most (but not all) classical 
bulges are quiescent and populate the red sequence of the di- 
agram. We presented evidence of different formation mecha- 
nisms for classical and pseudo-bulges. Classical bulges follow 
a correlation between Sersic index and bulge-to-total ratio, 
while pseudo-bulges do not. The latter are less concentrated 
than the former at fixed bulge-to-total ratio. The locus occu- 
pied by pseudo-bulges in the fundamental plane is different 
from that of classical bulges, and it is the same as that of 
discs. Furthermore, pseudo-bulges follow a mass-size relation 
different from that of classical bulges, and similar to that 
of bars. Altogether, this indicates that pseudo-bulges are 
formed through slow, non-violent processes, such as those 
expected from disc instabilities, whereas classical bulges are 
formed through violent processes, such as the merging of 
smaller units. We verified, however, a significant overlap in 
the properties of classical and pseudo-bulges, which sug- 
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Figure Al. Comparison between the results from test and final fits in twenty barred galaxies. Filled circles refer to tests where the 
input parameters were multiplied by two, while empty circles refer to tests where the input parameters were divided by two. Error bars 
are la uncertainties estimated by the fitting code. For B/T, typical error bars are shown, calculated through error propagation using 
the uncertainties in fie, and n. 



APPENDIX A: FIT SENSITIVITY TO INPUT 
PARAMETERS 

As mentioned in Sect. [3l image fitting often involves find- 
ing a global minimum in a very complex multidimen- 
sional parameter space. This means that, in principle (and 
depending on the algorithm used), fits can be very sensitive 
to the initial guesses given to the fitting code to start the 
first iteration. To avoid erroneous results, we have opted to 
produce and check each fit individually. Nevertheless, this 
alone does not mean that we find the most appropriate set 
of input parameters in all cases. 

In order to check how sensitive the results of our fits 
are to the set of initial guesses provided to budda, we have 
performed the following tests. We have randomly selected 
twenty barred galaxies, and produced test fits, in which the 
input values for /lo, h, fj,s, re and n used in the correspond- 
ing final fits are increased by a factor two (for ten of these 
galaxies) or decreased by a factor two (for the remaining 
ten galaxies). The choice for barred galaxies assures that 
the results from these tests concern typically difficult fits. 
Note also that a factor two is rather large, and that twenty 
galaxies correspond to 7 per cent of our barred galaxies. Fig- 
ure lAll shows the results from these test fits, as compared 
to those from the final fits. Error bars are la uncertainties 
estimated by budda. Its current version does not provide 
estimates of the uncertainty in B/T. Thus, for B/T, typi- 
cal error bars are shown, calculated through error propaga- 
tion using the uncertainties in fie, re and n, and assuming 
that the main source of error in B/T comes from determin- 
ing the total bulge luminosit y. Through similar tests with 
BUDDA, iDurbala et al.l (|2008l ) found comparable uncertain- 
ties in B/T (albeit somewhat smaller). 

One sees that, in most cases, the code converges to sim- 



ilar results, even for input parameters varying by a factor 
four. In fact, most results are consistent within the esti- 
mated uncertainties. All outliers come from the few fits that 
resulted in notably high values. Further more, these re- 
sults reinforce the findings in iGadottil (l2008l ). by indicating 
that disc parameters are particularly robust, and that B/T 
is a more robust parameter than n, which can have large un- 
certainties compared to its dynamical range. These results 
also argue in favour of our method to identify pseudo-bulges, 
at least for statistical studies such as ours, since He and re 
are also more stable than n. 

These results do not mean that any reasonable set of 
input parameters is appropriate. But they show that, by cre- 
ating input parameters through a careful inspection of each 
galaxy individually, we are able to achieve robust results. In 
fact, as noted above, wrong fits often result in significantly 
higher values (see leftmost panel in Fig. lAip . turning 
them generally easy to recognise in our procedure. Only 2.5 
per cent of our final fits have x^ > 2. 



